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1 Introduction 



Supersymmetric Yang-Mills quantum mechanics (SYMQM) shares many prop- 
erties of advanced field theories thereby it provides a useful laboratory to 
study explicitly some of their properties (HE]- It emerges from the dimen- 
sional reduction (in space) of the full supersymmetric Yang-Mills field the- 
ory defined in the D dimensional space-time. Depending on D it covers a 
wide range of interesting phenomena. D = 2 system is exactly soluble with 
manifestly supersymmetric continuous spectrum and analytically calculable 
Witten index QjJ |3] . For D = 4 the model is nontrivial and posesses both 
localized and non localized eigenstates, which have been constructed only re- 
cently in all fermionic sectors. On the other hand in the zero fermion sector 
it reduces to the zero volume limit of the pure Yang-Mills theory well studied 
in the context of lattice field theory jHElE]- Finally the D = 10, SU(oo) sys- 
tem, with its threshold bound state and the continuum of scattering states, 
is considered as a model of M-theory [7j and has attracted a lot of interest 
(for recent review see e.g.jH]). 

In this talk I would like to review the current status of the programme 
which attempts to solve these models in various dimensions and for various 
gauge groups 0. The four dimensional (D = 4) system with the SU(2) gauge 
group, where the main progress has been achieved, will be discussed in the 
next Section. In Section 3, a new approach to the D = 10, SU(2) model, 
together with some preliminary results for the purely bosonic sector will be 
presented. 

2 D=4 Supersymmetric Yang-Mills quantum 
mechanics 

Reduction of the D dimensional supersymmetric SU(2) Yang-Mills field the- 
ory to a single point in the d — D — 1 dimensional space leads to the quantum 
mechanical system which for D = 4 is described by nine bosonic coordinates 
x l a (t), i — 1, 2, 3; a = 1, 2, 3 and six independent fermionic coordinates con- 
tained in the Majorana spinor i[>%(t), a — 1, ...,4. Equivalently (in D=4) one 
could impose the Weyl condition and work with Weyl spinors. Hamiltonian 
reads j2] 

H = H B + H F , 
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H B = -PaPl + —t abc t ade X l b X 3 c X l d X 3 el (1) 

H F = % j-e abc ^T k ^ b x k , 

where ip T is the transpose of the real Majorana spinor, and V in D=4 are 
just the standard Dirac a matrices. 

After the reduction of the three dimensional space to a single point, the 
rotational symmetry of the original theory becomes the internal spin(3) sym- 
metry. It is generated by the angular momentum 

T = e^ k (x{p k a - i^E»"ty B ) , (2) 

with 

5? fc = -i[T i ,r*]. (3) 
The model has also a gauge invariance with the generators 

Ga = e abc [x k b p k - ^6 T ^c) , (4) 
and is invariant under the supersymmetry transformations generated by 

Qa = (rVa)aPa + igtabc^^a) (5) 

The bosonic potential (written in the vector notation in the color space) 

V = £v jk @ x x k )\ (6) 

exhibits the famous flat directions responsible for a rich structure of the 
spectrum. 

2.1 Original approach 

We shall now review a simple way to calculate automatically algebraic ex- 
pressions for matrix elements of a wide class of quantum observables and 
to construct numerically the complete spectrum of a given system j^]. The 
method implements literarily in the computer the rules which govern the 
quantum world. 
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Quantum mechanics inside a PC 

For a moment consider a single bosonic degree of freedom. Generalization 
to fermions and to more variables readily follows. Action of any polynomial 
observable can be easily implemented in an algebraic program if we use the 
discrete eigenbasis of the occupation number operator a^a 

{\n>}, \ n >=^=^r\0>. (7) 



n 

Since the bosonic coordinate and momentum operators are 

x = —=(a + a^), p— — =(a — ah, (8) 

y/2 y i^r ' v ; 

a typical quantum observable can be represented as the multiple actions of 
the basic creation and annihilation operators 1 . 

A quantum state is a superposition of arbitrary number, n s , of elementary 

states |n > 

|st >= E^a/|n (/) >, (9) 
and will be represented SIS db Mathematica list 

st = {n s , K ...,a ns }, {n^}, . . . , {n<->}}, (10) 

with n s + 2 elements. The first element specifies the number of elementary 
states entering the linear combination, Eq.fjUJ). The second element is the list 
itself and contains all complex amplitudes aj,I = 1, . . . ,n s . Remaining n s 
sublists give the occupation numbers of elementary, basis states. According 
to this convention an elementary state |n > is represented by {1, {1}, {n}}. 

Next, we implement basic operations defined on states: addition, multi- 
plication by a number and the scalar product. They are simply programmed 
as definite operations on Mathematica lists transforming them in accord with 
the principles of quantum mechanics. Creation and annihilation operators 
are then defined as a list- valued functions on above lists. According to Eq. (JH]) 
the action of the position and momentum operators becomes also defined. 
Then we define any quantum observable: hamiltonian, angular momentum, 
generators of gauge transformations, supersymmetry generators, etc. 

Further procedure is now clear: given a particular system, define the list 
corresponding to the empty state, then generate a finite basis of N cut vectors 

1 The method can be also extended to non polynomial potentials. 
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and calculate matrix representations of the hamiltonian and other operators 
using above rules. Next, the complete spectrum and its various symmetry 
properties is obtained by the numerical diagonalization. 

By studying the cutoff dependence of the spectrum one can estimate the 
systematic errors induced by restricting the Hilbert space. In many systems 
studied so far convergent results were obtained before the basis grew too 
large. 

Creation and annihilation operators for SYMQM 

Super symmetric hamiltonian is polynomial in momenta and coordinates, 
hence above idea can be readily applied. To this end rewrite bosonic and 
fermionic variables in terms of the creation and annihilation operators of 
simple, normalized harmonic oscillators 

=*2. a p x f } = c p,*=i,2, (ii) 

such that the canonical ( ant i) commutation relations 

[xlp k b ]=iS ik S ab , {C,<} = ^f- (12) 
are preserved. Standard extensions of Eq. (JSJ) for bosonic variables read 

< = ^K + a%), pl = -^K-a^). (13) 

For fermions, the following representation for a quantum hermitean Majorana 
spinor was used 



/ 1 + i 



^ifa ~ fa ~ fa if a 
~f a if a if a ~~ fa 

K-ifl-fl + ft+iff) 



(14) 



The basis and the cutoff 

The complete Hilbert space is spanned by all independent polynomials of 
creation operators a\ and f^ acting on the empty state 

|(0, 0, 0), (0, 0, 0), (0, 0, 0), (0, 0, 0), (0, 0, 0) >, (15) 

which in the Mathematica "representation" reads 

{1, {1}, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}}}. (16) 



5 



By construction, the first three vectors (in color) specify bosonic, and the 
last two fermionic, occupation numbers. In practical applications we shall 
work in the restricted Hilbert space containing at most B bosonic quanta in 
total. Hence the gauge and rotationally invariant cutoff N cut is defined as 

£i,&<4a|, = B < B max = Ncut. (17) 

Since the Pauli principle admits only six Majorana fermions in this system, 
there is no need to restrict the fermion number F. 

Local gauge invariance is taken into account by constructing only the 
physical, i.e. gauge invariant basis. To create all independent, gauge invari- 
ant states at fixed F and B consider all possible contractions of color indices 
in a creator of (F, B) order 

^■■WX* ( 18 ) 

for all values of the spatial indices i and o. All color contractions fall naturally 
into different gauge invariant classes. Creators from different classes differ 
by color contractions between bosonic and fermionic operators. For example 

WWW (is) 

and 

are in different gauge invariant classes. Creators of odd order are constructed 
with one triple contraction. For example 

e^afaKaf/jVf (21) 

and 

e cde a^afa^afrjf^ (22) 

also belong to different gauge invariant classes. To select linearly independent 
states we used again the rules of "quantum algebra". Therefore, at fixed 
F and B the final procedure is as follows: (i) identify all gauge invariant 
classes of creators, (ii) loop over all values of spatial indices and for each 
ii, ib, 0i, gf create corresponding state from the empty state, Eq.(|I5j). 
(Hi) identify and reject linearly dependent vectors, (iv) orthonormalize the 
remaining set of states. 
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Figure 1: Spectrum of D=4 supersymmetric Yang-Mills quantum mechanics. 

Given the basis it was then a simple matter to calculate automatically 
matrix representation of the hamiltonian and other observables. 

Early results 

In Fig. 1 we quote the spectrum obtained in Ref . |H] in all four independent 
fermionic sectors 2 . The cutoff in this calculation was varying between B max = 
8 for F = and B max = 5 for the largest (i.e. most difficult) F = 3 sector. 
The angular momentum of a sample of states is also shown. It was deduced 
from the degeneracy of the SO (3) multiplets and independently from the 
explicit representation of the angular momentum in our basis. 

Later on this results were upgraded, with the considerable numerical ef- 
fort, to B max = 8 in all sectors Fig. 2 is the outcome of the latter cal- 
culations which took a couple of months of a 600 MHz ALPHA workstation. 
The cutoff dependence displayed in this Figure confirmed and quantified the 
general expectation that the system has both continuous and discrete spec- 
trum. Discrete, localized states (which manifest themselves here as quickly 
convergent levels) appear for F = and F = 1 while the non localized states, 

2 The hamiltonian has also the particle-hole symmetry. 
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Figure 2: Cutoff dependence of the first few eigenenergies of the hamiltonian 



with substantially slower convergence, are seen in the "fermion rich" sectors 
with F = 2 and 3. This was expected from the flat valley nature of the 
potential, Eq.lJHjl. combined with the supersymmetry [TTJ. A series of other 
results about the supersymmetric structure of the system, Witten index, etc. 
was also obtained. 

Further progress was achieved recently by exploiting analytically the sym- 
metries of the system. 



2.2 New developments 

Separation of variables 

Above results have been beautifully confirmed and extended with the aid of 
the "almost analytical" approach pioneered by Savvidy [12] and developed by 
van Baal in slightly different context ^3] • Decomposing solution of the nine 
dimensional Schrodinger equation, in the F = 0, J = and F = 2, J = 
channels into covariant tensors, he reduced the problem to a numerically 
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tractable set of coupled ordinary differential equations. When adapted to 
our case his method can push the cutoff as high as B max = 39 in these two 
channels, see Fig.3. Now the discrete, localized, and quickly convergent with 
the cutoff states with F = are clearly seen, their energy determined with 
a very high precision. Moreover, the intricate nature of the solutions with 
two fermions is also evident. The flat lines signal again the localized bound 
states of two gluinos, while slowly falling with the cutoff levels correspond to 
the non localized states from the continuous spectrum. Some of the bound 
states in F = and F = 2 sectors have degenerate energies as required 
by supersymmetry. Other supersymmetric partners must then be located in 
fully fermionic (F = 1,3) sectors for which the method has not been yet 
generalized. 



F=0 J=0 F=2 J=0 

20 40 20 40 




-0.5 



20 40 20 40 

B 

Figure 3: High cutoff results from van Baal approach. 
High cutoff solutions in all channels 

Similar high precision results are now available for all angular momenta and 
in all fermionic sectors [T3|. To this end we have extended the recursive 
approach first applied to the D = 2 system [3]. In the two dimensional case 
recursions, which relate the matrix elements in a bigger basis with those 
in the smaller basis, were simple enough to be solved analytically. For the 
D = 4 this is not the case any more. However they speed up dramatically 
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computations allowing to reach high precision spectra in all fermionic sectors 
and all angular momentum channels. Second, rotational invariance is now 
fully employed, by projecting into channels of given angular momentum. 
Consequently there is no need to diagonalize very large matrices. Both these 
improvements result not only in the more precise spectra, but in much more 
satisfactory understanding of the supersymmetric structure of the D = 4 
system, including complete classification of all supermultiplets. 
Analytical study of continuous spectra in cut systems 
Finally, in a recent development a new scaling in the continuous spectra of 
a cut quantum systems was found jTSj. This provides an important tool for 
studying the continuum limit of non localized states and their interactions. 

To summarize this section, the technology to study quantitatively D = 4, 
SU(2), SYMOM is now available and one can readily attack more subtle 
problems like scattering ^H] or thermodynamics of this system ^H] • 

We now turn to, much more difficult, D — 10 model whose quantitative 
study with presently discussed methods is just beginning (see [21112] for other 
attempts). 

3 TWS scheme 

In higher dimensions only the original approach is available and the spectrum 
of lowest zero- volume glueballs was obtained for all 4 < D < 10 in the F = 
sector ^7j. The number of states grows rapidly with D and consequently at 
D = 10 only B < 4 were reached in practice. In fermionic sectors situation 
is yet worse. The difficulties are two-fold: a) bases become too large to 
store, and b) time required for calculating the full hamiltonian matrix grows 
prohibitively large. The TWS scheme, presented below, offers substantial 
reduction of both parameters. 

Consider the (B,F) sector of the Hilbert space, i.e. the one with B bosons 
and F fermions. The basis vectors before orthonormalization can be uniquely 
characterized by the compound index 

/={<7,i\...,iVV--,^}, (23) 

since given / the actual state can be readily created by the compound creator 
from the gauge invariant class C (c.f. nT?ll2*2"j) . Many properties of the mixing 
matrix and of the hamiltonian matrix can be deduced from the index / 
alone without need to create the actual state. Therefore we shall abandon 
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where possible the lengthy Mathematica representation of states and use 
instead much more compact indices (|23|) . This virtually eliminates the first 
problem. Second difficulty is also reduced by observing that for polynomial 
hamiltonians most of the matrix elements are zero and many of non-zero 
ones are actually the same. We shall therefore expose the sparse nature of 
the hamiltonian matrix and identify equal matrix elements. This can be done 
using indices (J2HJ) only, and will lead to the new classification of basis states. 

We begin with the properties of the mixing matrix Mjk =< I\K >, also 
referred to as the norm matrix. 

3.1 Norm matrix and new classification of states 

For the sake of simplicity we restrict the analysis to the F = sector. Gener- 
alization for fermions is possible. Consider the generic element of the mixing 
matrix 

M IK =< I\K >=< C I ,i 1 ,...,i Bl \C K ,k\...,k Bli > . (24) 

Obviously only states which contain the same quanta can mix, therefore 

M IK ^0 if B I = B K = B,{i l ,...,i B }~{k\...,k B }, (25) 

where the " ~" means that the values of all { i } and { k } indices coincide up 
to a permutation. On the other hand states from different gauge invariant 
classes can mix in general. 
Words 

This simple observation suggests to group all states with given B into words. 
A word W is defined as a set of all states which contain the same quanta. 
For example in a B = 4 sector states {0,1,2,2,5} and {D, 1,2,5,2} belong to 
the same word while states {C',4,3,3,8} and {D', 1,2,3,5} are from different 
words. All gauge invariant classes, which can be constructed for particular 
set of spatial indices, are included in a given word by definition. Hence a 
word is uniquely labelled by an ordered set of B integers, e.g. W=(1225) 
defines the word considered above. How many states are in a word? As 
many as there are inequivalent permutations of spatial indices which give 
linearly independent states in all contributing gauge invariant classes. Now 
Eq.fl23j) can be stated as the following 
THOREM 1: Only states in the same word mix. 

Hence the whole mixing matrix splits into blocks of small (Nw x Nw) 
mixing matrices for each word W with Nw being the number of states in a 
word W . 
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Types 

Words come in different types. A type T characterizes the pattern of indices 
irrespectively of their actual value. For example words (1225) and (3448) 
belong to the same type which we denote as [211], and (1234) or (2227) 
belong to the two different types [1111] and [31] respectively The number 
of words in a given type, is equal to the number of realizations of the 
pattern by different values of the indices. For our first example (d — D — 1) 

jV 211 = d(d-l)(d-2)/2, (26) 

and similarly for other types. Types can be identified with the partitions of 
B. Therefore the total number of types in a given bosonic sector is just equal 
to the number of partitions, Pb, of B. 

Since neither the number of states in a word nor the number of types in a 
sector depend on the dimensionality of the problem, the main complexity lies 
in the number of words in a given type, e.g. Eq. (J26j) which grows rapidly with 
d for larger B. However the big simplification occurs due to the following 
consequence of the Wick theorem. 

THEOREM II: All mixing matrices in the same type are equal. 
To see this rewrite Mjx as the sum over all contractions 

M IK =< I\K >=< 0\a h . . . a iB a ] kl . . . 4j0 >= E COTfrac a orM ILJ i . fc ., (27) 

where the product on the right hand side denotes symbolically given contrac- 
tion. It is evident, even from this schematic expression, that Mik does not 
depend on the particular values of the indices, but depends only on their 
configurations or patterns. Hence the theorem follows. 

Therefore the calculation of the full mixing matrix is simply reduced to 
calculating Pg small matrices and using the same copy for all words in a given 
type. Moreover, the tedious orthonormalizaton procedure also decouples into 
small blocks, corresponding to words, and has to be done only for one word in 
a type with the rest being the exact copies of the results from the first word. 
This simplification prompts us to introduce the new scheme of organizing 
states in a basis - the TWS scheme which groups states according to the 
following hierarchy: 

S {totes) — > W{ords) — > T{ypes) — > Sectors. (28) 

With this organization any state in a basis is labelled by four integers 

| J >= \b,t,w,s >, (29) 



12 



AB 


rule 


action 






1 


diagonal 







2 


2^2 new 






3 


2 ^2 old 






4 


add 2 new 




2 


5 


add 2 old 






6 


add 2x2 new — 


new 


4 


7 


add 2x2 new - 


- old 




8 


add 2x2 old — 


old 



Table 1: Effective transitions induced by the bosonic hamiltonian (P) in nine 
space dimensions. 

with b, t, w, s enumerating sectors, types, words and states respectively. 

It turns out that the structure of the hamiltonian is also much more 
transparent in this representation as discussed below. 

3.2 SYMQM hamiltonian 

As above we shall consider only bosonic sectors (B, F) = (B, 0) where only 
Hb of Eq.((H) contributes. Since Hb contains terms of the second and fourth 
order in bosonic creation and annihilation operators, it can only induce tran- 
sitions with AB = 0, 2, 4. Each transition can be interpreted as the effective 
operator creating/annihilating 0,2 or 4 indices in our TWS representation. 
They are listed in Table 1 and will be referred to as rules since each transition 
implies a rule how to generate from the initial type/word the final type/word 
with nonzero matrix element. For example AB = transition, labelled as 
rule 2, annihilates two indices in the initial state and creates 2 indices with 
a common value which is new, i.e. does not exist in the initial state. Rule 3 
has similar action but the common value of the created pair coincides with 
one of the values already existing in the initial state. Rule 7 adds two pairs of 
identical indices one being new and one already existing, etc. It follows from 
Eq.© that we cannot create/annihilate a pair of different indices. Similarly 
there are no terms creating four the same indices. 

We emphasize that above classification should be used only to identify 
all nonzero matrix elements of the hamiltonian. The actual values of matrix 
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4 5 
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5 
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1 2 

1 3 2 

2 123 2 
3 1 2 

1 3 2 

2 2 3 123 2 
2 12 

13 2 

2 2 12 

2 12 

1 



Table 2: Structure of the bosonic hamiltonian in the TWS representation, b 
and t are defined in Table 3. 



elements will be calculated within the original approach. 

Table 2 shows the resulting structure of the hamiltonian in the TWS 
representation for the first four sectors. Labelling of rows and columns is ex- 
plained in Table 3. The hamiltonian is block tri-diagonal in bosonic sectors 
(odd and even sectors are decoupled bacuse of the parity conservation). Each 
nontrivial sector-to-sector block is again sparse with only few nonzero sub- 
blocks connecting initial and final types according to the rules listed in Table 
1. For simplicity they are labelled by these rules. Each such type-to-type 
subblock is again a N?f x N-n matrix with entries connecting different words 
of the initial and final types. Again this matrix is in general sparse - from 
given initial word one can reach, with a given rule, only few final words. This 
is not the end, the entries of this matrix are again matrices labelled by the 
individual states belonging to initial and final words. Only at this level the 
orthogonalization of basis states enters. Since elementary states in a single 
word are not orthogonal, the word-to-word blocks must be corrected by the 
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Table 3: Labelling of sectors and types in Table 2. 
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406 


666 
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5 


6 


36 


126 


336 


756 


1512 


2772 


6 


56 


220 


680 


1771 


4060 


8436 


16215 


7 


21 


180 


855 


2976 


8478 


20952 




8 


126 


714 


3045 


10521 









Table 4: Sizes of bases in the F = sector for space-time dimensions 4 < 
D < 10. N s is the number of basis vectors with given number of bosonic 
quanta, B. 

square roots of the corresponding norm matrices. This is, however, much less 
time consuming than the global orthonormalization. 

Moreover one can show that the word-to-word subblocks in a given type- 
to-type block are independent of the initial word. Similarly to the mixing 
matrix case, it suffices to calculate transitions from one initial word in each 
type and substitute them for other initial words. The proof is more technical 
and will not be given here. 

The whole structure, even though little involved at first glance, can be 
readily formalized and leads to much faster calculations of the whole matrix. 

3.3 Preliminary results 

The complexity of the problem is summarized in Table 4 where dimensions 
of the Hilbert space for various D and B are displayed. All bases were 
generated and orthogonalized with the original approach (c.f. Section 4.1). 
Mathematica programs took about 2 days to produce the last entry, thereby 
the cutoff 5 = 7 — 8 seems to be a practical maximum for the brute force 
method at D = 10. Calculation and diagonalization of the hamiltonian is yet 
more demanding - for D=10 one can reach only B = 4 within the reasonable 
computing time [T7] . 
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The TWS approach allows to reach one to two orders of magnitude bigger 
bases, which at D=10 translates for B max = 6 — 8. This rather modest, in 
terms of B, improvement is nevertheless quite relevant for the precision of 
the energy levels. Fig. 4 shows the first five eigenenergies as a function of 
the cutoff. Solid lines are drawn just to guide the eye. Levels are labelled by 
corresponding SO (9) representations. Since our cutoff respects the rotational 
invariance, resulting spectrum of a cut hamiltonian should have full SO (9) 
symmetry 3 , and indeed we observe first few SO(9) multiplets as indicated 
in the Figure 4 . The ordering of the levels (singlet, tensor, singlet) turns 
out to be the same as in the well known D = 4 case [S] and in fact remains 
unchanged for all 4 < D < 10 ^7j . The energy of the lowest singlet state 
has already converged within 3%. For second and third state the relative 
change with B dropped from 13% to 7% and from 17% to 10% respectively. 
Cutoff B max = 8 can be reached with the reasonable computing effort. 

120 
100 
80 
E 60 
40 
20 



2 3 4 5 6 7 

B 

Figure 4: Lower levels of the spectrum of the pure Yang-Mills D=10 system. 



3 Only in the purely bosonic sector. 

incidentally this provides one of the tests of the whole TWS scheme. 
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4 Summary and outlook 



In the last year there has been a substantial progress in solving supersym- 
metric Yang-Mills quantum mechanics with the SU(2) gauge group. In four 
dimensions we can reach such high cutoffs that the restriction of the Hilbert 
space becomes irrelevant for most purposes. Precise spectrum in all channels 
has been obtained and the intricate pattern of continuous and discrete states 
found and clarified. Supersymmetric structure of the system can be now 
studied in detail. In particular the complete classification of supersymmetric 
multiplets is now available. 

In the ten dimensional system developments are obviously slower but nev- 
erheless steady. In the bosonic sector one can now reach the cutoffs (hence 
the precision) comparable to the D = 4 model studied with the original tech- 
nique. The new TSW scheme is yet to be generalized to fermionic sectors. 
The ongoing puzzle of nonconserved massless Major ana fermions calls for 
further study. On the other hand reported here progress in the four dimen- 
sional case opens some new possibilities also for the ten dimensional system. 

Acknowledgments I thank the Theory Group of the Max-Planck-Institut in 
Munchen for their hospitality. This work is supported by the Polish Com- 
mittee for Scientific Research under the grant no. PB 2P03B09622, during 
2002 - 2004. 

References 

[1] M. Claudson and M. B. Halpern, Nucl. Phys. B250 (1985) 689. 

[2] M. B. Halpern a nd C. Schwartz, Int. J. Mod. Phys. A13 (1998) 4367, 
|hep-th/9712133 

[3] M. Campostrin i and J. Wosiek, Phys. Lett. B550(2002) 121, 
|hep-th/0209i"40l 

[4] M. Liischer, Nucl. Phys. B219 (1983) 233. 

[5] M. Liischer and G. Minister, Nucl. Phys. B232 (1984) 445. 

[6] P. van Baal, Acta Phys. Polon. B20 (1989) 295. 



18 



[7] T. Banks, W. Fishier, S. Shenker and L. Susskind, Phys.Rev. D55 
(1997) 6189; |hep-th/9610043| 



[8] W. Taylor, Rev. Mod. Phys. 73 (2001) 419, |hep-th/0101126 



[9] J. Wosiek, Nucl. Phys. B 644 (2002) 85, |hep-th/0203TT6 

[10] J. Kotanski and J. Wosiek, Nucl. Phys. 5(Proc. Suppl.)119 (2003) 932, 
|hep-lat/0208067| 

[11] B. de Wit, M. Liischer and H. Nicolai, Nucl. Phys. B320 (1989) 135. 
[12] G. K. Savvidy, Phys. Lett. B 159 (1985) 325. 

[13] P. van Baal, The Witten Index Beyond the Adiabatic Approximation, 
in: Michael Marinov Memorial Volume, Multiple facets of quantization 
and supersymmetry, eds. M. Olshanetsky and A. Vainshtein (World Sci- 
entific, Singapore, 2002), pp.556-584, hep-th/0112072. 



[14] M.Campostrini and J. Wosiek, in preparation 

[15] M. Trzetrzelewski and J. Wosiek, Jagellonian University preprint, 
|hep-th/0308007| 

[16] U. H. Danielsson, G. Ferretti and B. Sundborg, Int. J. Mod. Phys. All 
(1996) 5463, |hep-th/960308H 

[17] J. Wosiek, Supersymmetric Yang-Mills quantum mechanics, in: Pro- 
ceedings of the NATO Advanced Research Workshop on Confinement, 
Topology and Other Non-Perturbative Aspects of QCD, eds. J. Green- 
site and S. Olejnik, Kluwer AP, Dordrecht, 2002, |hep-th/0204243| 

[18] K. Becker and M. Becker, Nucl.Phys.B506 (1997) 48-60, 



hep-th/9705091 



[19] D. Kabat, G. Lifschytz and D. A. Lowe, Phys. Rev. D64 (2001) 124015, 
|hep-th/010517l| 



19 



